Decoherence and Quantum Walks: anomalous diffusion and ballistic tails 
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The common perception is that strong coupling to the environment will always render the evo- 
lution of the system density matrix quasi-classical (in fact, diffusive) in the long time limit. We 
present here a counter-example, in which a particle makes quantum transitions between the sites 
of a d-dimensional hypercubic lattice whilst strongly coupled to a bath of two-level systems which 
'record' the transitions. The long-time evolution of an initial wave packet is found to be most 
unusual: the mean square displacement n 2 of the particle density matrix shows long-range ballitic 
behaviour, with (n 2 ) ~ t 2 , but simultaneously a kind of weakly-localised behaviour near the origin. 
This result may have important implications for the design of quantum computing algorithms, since 
it describes a class of quantum walks. 

PACS numbers: 05.40.Fb, 03.65.Yz, 03.67.-a 



One can think of the trajectory of a quantum parti- 
cle hopping between 2 nodes A and B on some lattice 
or 'graph', as a 'quantum walk', in which the amplitude 
to go from A to B is given by summing over all possi- 
ble paths (or 'walks') between them. Amusingly, such 
walks can also describe the time evolution of quantum 
algorithms, including the Grover search algorithm and 
Shor's algorithm. One can find explicit mappings be- 
tween the Hamiltonian of a quantum computer built from 
spin- 1/2 'qubits' and gates, and that for a quantum par- 
ticle moving on some graphpJQ- Each graph node repre- 
sents a state in the system Hilbert space, and the system 
then walks in 'information space'. This mapping is most 
transparent for spatial search algorithms with the local 
structure of the database. Amongst the graphs so far 
studied are 'decision trees ' d H and hypercubesd @ ; 
quantum walks on other graphs, and their connection to 
algorithms, were recently reviewed 0. 

The quantum dynamics between two sites A and B on 
a given graph is often much faster (sometimes exponen- 
tially faster) than for a classical walk on the same graph 
[J • It nas been argued that quantum walks may 
generate new kinds of quantum algorithm, which have 
proved very hard to find. Several recent papers have 
also considered experimental implementations of quan- 
tum walks for quantum information processing^, Il0| : 
some involve walks in real space, whereas others are 
purely computational (eg ., a walk in the Hilbert space 
of a quantum register[l(J). Many experiments over the 
years, particularly in solid-state physics, have also been 
implicitly testing features of quantum walks. 

As always, the main problem confronting any quantum 
algorithm is environmental decoherence - the gradual en- 
tanglement of the system with the 'environment' means 
that phase interference effects are gradually lost, in mea- 
surements performed on the system alone. It is generally 
assumed that the system dynamics will then show classi- 
cal diffusion at long time s \T\\ , at least if the environment 
is at or near equilibrium 12]. This 'folk theorem' is sup- 



ported by results on many modelsll3 | (e xcept for certain 
very unusual 1-dimensional systems |l4|). Recent inves- 
tigations of decoherence effects on quantum walks |l Si- 
ll!! give similar results, although in one investigation 
of random walks driven by coin-tosses [17], non-classical 
behavior was found. In these recent investigations, the 
decoherence mechanism was either (i) an external noise 
source (ii) a coupling to a set of tossing 'coins'; or (iii) a 
coupling of the coins to a heat bath. In solid-state and 
atomic qubits systems, the heat bath is modelled either 
by a set of oscillators (representing delocalised modes 
like phonons, photons, or electrons), or by a set of '2- 
level systems', or 'TLS' (representing localised modes like 
defects, topological disorder, or nuclear and paramag- 
netic spins). Both are important in experiment; TLS are 
particularly important for decoherence in mag netic[2(|, 
superconducting[2l|, and conducting 22] qubit systems, 
and tend to dominate at low temperature. 

In this paper we consider a class of quantum walk mod- 
els having a very unusual dynamics- not only is the long- 
time behaviour not classically diffusive, but a part of the 
single-particle reduced density matrix always continues 
to show coherent dynamics. These models are very rele- 
vant to solid-state quantum information processing sys- 
tems, since they involve a TLS bath- we couple a quan- 
tum particle moving on a graph to a bath described by 
a set {<Tfe} of TLS, written as spin-1/2 Pauli spins (with 
k = 1,2, ...N). We first describe the dynamics of these 
models, and then their physical interpretation. 

Quantum Walker: For definiteness we choose a d- 
dimensional hypercubic graph for the walking particle 
(our main conclusions do not depend on this assump- 
tion), with the 'bare' Hamiltonian 



Ho = A a J2 (4cj + h-c.) 



(1) 



Here cj creates the particle on site i, and < ij > restricts 
the dynamics to nearest neighbor hopping. The particle 
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moves in a Bloch band with dispersion relation e e (k) = 
2A ^2 , cos(k l _ t a ) and bandwidth W = AdA . Here 
a is the lattice constant, and k the d-dimensional mo- 
mentum. Henceforth we measure all distances in units of 
a Q , and label lattice sites by a lattice vector n. 

For this quantum walker, the solution of Schrodinger's 
equation is standard. Thus a particle initially localized 
at the origin, with wave-function ip n (t = 0) = 5 n0 at 
t = 0, evolves to ip n (t) = L~ d £ k e '[k n- £o (k)]t at a later 
time, where L is the linear system size. The probability 
distribution P° (t) — \ip n (t)\ 2 is then 



(2) 



where J n (z) is the n-th order Bessel function. The 
continuum-space limit is recovered by considering a 
broad Gaussian initial wave-packet, initially centred at 
the origin, of form ip n (t = 0) « (l/^R) d/2 e -™ 2/2i?2 
with R 3> 1. Then for later times 



R 2 



7r(i? 4 + z 2 ) 



d/2 



-n 2 R 2 /(R i +z 2 ) 



(3) 



As expected, a purely quantum-mechanical evolution 
gives P§ {t) oc l/t d and (n 2 ) oc i 2 at long times. 

Environmental decoherence: Coupling the quan- 
tum walker to an environment is supposed to change the 
long-time evolution to classical diffusion, characterized 
at long times by Pffl(t) oc l/t d / 2 and (n 2 cl ) oc t. We 
certainly expect this for models in which the particle co- 
ordinate is coupled to an Ohmic oscillator bath, but we 
now examine the effect of a coupling between the particle 
and a TLS bath. On its own, this bath has a Hamiltonian 

Ktls = Efc h fc ' o-fc + J2kk' Kk'VkVk-i where the { h fe} 
are fields acting on each TLS, and the V£jS describe in- 
teractions between them. Typically the are very 
small, and lead only to a very slow dynamics of the TLS 
bath, so we shall drop thempjj. Various couplings of the 
bath to the walker are possible, but we are interested in 
those which monitor transitions of the walker, ie., those 
triggered when the particle hops between nodes. We can 
then distinguish 2 important limiting cases: 

(i) the TLS bath is acted on by only weak external 
fields, which we then neglect. Now assume that each time 
the quantum walker hops it can flip the fc-th TLS <Jk with 
amplitude ctk- We can write the effective Hamiltonian as 



(ii) The TLS bath is polarised by strong external field 
hfc = h. Defining the unit vector {z} along the axis of 
this field, and the total polarisation M = ^2 k af. of the 
TLS bath with respect to this axis, we see that in this 
strong field limit, only bath transitions which conserve M 
are allowed. In this case one has an effective Hamiltonian 

Km = A £ [c\c [P_ M e*£* ak ° l P M \ + H.c] (5) 

where Pm projects the state of the TLS bath onto the 
subspace with polarization M. We have dropped the 
large Zeeman term ^ fc h-erj, from this Hamiltonian, since 
it is now just an M-dependent constant. 

We now proceed with the solution for the probability 
distribution P n o(t). We shall look in detail at the first 
model Q above, and then comment on the second one. 
The form of 10} is a simple generalisation of a Hamilto- 
nian Tt = A a {t x cos[J2k a kVk] + H.c.}, which describes 
one limiting case of the interaction of a single qubit r 
with a spin bath. The density matrix of this model is 
given exactly as aphase average over the propagator of 
the 'bare' qubit jU 0]> an d one can use precisely the 
same technique to write the solution for 10} . Thus, for 
the initially localised state tp n (t = 0) = <5 n o, assuming 
the strong decoherence limit described above, one finds 
the solution at time t as 



Pno(t) 



2,T 



2^ 



II J n»( ZC0Si P) 



(6) 



and similarly for the initially broad wave-packet one gets 



Pn (t) 



dip R d e - n R /( R +z cos ^ 



o 27r [ir(R 4 + z 2 cos 2 ip)} 



d/2 



(7) 



We will rederive this result using a rather different 
method at the end of the paper. Given the strong cou- 
pling to the TLS environment, one intuitively expects 
classical diffusive dynamics at long times. Surprisingly, 
the actual evolution is radically different. Consider first 
the probability at long times of finding a particle back 
at the origin, Pqo(z — > oo), in Eq. ©. The asymp- 
totic expansion for the Bessel function Jo(zcosco) ~ 
\plpK cos(2 cos if — it /4) / (z cos ip) is not possible because 
cos ip — > for if — > ±7r/2. In fact, in the t — > oo limit the 
dominant contribution (for d > 1) comes from tp w ±vr/2. 
Then 



H 




H.c. 



(4) 



In what follows we will assume that the individual 
are small but that the number TV of TLS is so large that 
k = J2k a k 1> i- e - hopping events are accomplished 
by simultaneous transitions in a large number of TLS. In 
other words, we look at the case of strong decoherence. 



Poo(z 



oo) ~ — 

7T 



d V 4 d {z V ) = A , (8) 



where Ad = (27r) — 1 J dxJQ d (x) is a constant (in d = 1 
there is an additional ln(2A c t) factor). This result is 
already rather peculiar since in d > 2 the decay of Pqo (f) 



x 



t -d/2 } 

and in the ideal, or ballistic, quantum propagation, P° m oc 



is integrable both in the classical diffusion, P ( 



(cZ) 
00 
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t . We get qualitatively similar answers for the broad 
initial state Q, where 

A K>2 — d 

P 00 (z^cx,)«^_, (9) 

and A d = 2ir- ( - 1+d / 2 '> J °° dx/(l+x 2 ) d / 2 . The dependence 
on the initial wave-packet spread in the long time-limit 
is another unusual feature of the solution. 

From the divergence of the total time spent at the ori- 
gin t = /q - * 00 (1tPq{t) oc A~ 1 ln(A G i) — > oo one's first 
suspicion is that the strong environmental decoherence 
is causing some kind of quasi-localization of the particle, 
analogous to weak localisation in solid-state physics. It 
then comes as an astonishing paradox that a calculation 
of the mean-square displacement from JfjJ) gives 

< ((n(t) - n(0)) 2 >=12^n 2 P n0 (z) = ^(A o i) 2 , 

n 

(10) 

which is only a factor of two smaller then the coher- 
ent quantum evolution! Thus the solution shows quasi- 
localsation near the origin, coexisting with coherent bal- 
listic dynamics at large distances. 

Having both Poo(i) cc 1/* an< ^ E n n2 ^oW oc t 2 at 
the same time is obviously inconsistent with the sim- 
ple scaling form t~ d f(n 2 /t 2 ). The solution to the para- 
dox requires a more complex shape for the distribution 
function, which we show in Fig. 1 and derive here for 
the Gaussian initial state. We introduce new variables 
r = n/R and u = z/R 2 to simplify the integral in to: 

f I \ d/2 f 2 ^ dip e -r 2 /(i+« 2 cos 2 v ) 
P ro[u) = -— j / — — — , 1d/2 • (11) 

It is straightforward at this point, by considering the 
long-time limit u > 1, to derive the following relations 
for the intermediate 

^^^S^A^T' (R«n«A ot /R), 

(12) 

and large (n — > oo) length scales 

P „W rj \ ( R \ 1 c -n 2 (R/2A n t) 2 (13) 

As expected from (0 , one has an increased probability of 
finding a walker at the origin. Since the power-law decay 
l/n d ^ 1 is not integrable, the normalization integral and 
(n 2 ) are still determined by the parameter z, but the 
probability of being at the origin is enhanced over that 
at a distance A t/R by a factor (A i/i?) d_1 . 

Not surprisingly, the thermodynamics of this sys- 
tem is also peculiar. The partition function Z = 
f 2v (dip/2Tr)lQ (2A Q cos ip/T) is leading to a free energy 
in the low-temperature limit T <gC A Q given by F(T) ~ 
-T(d/2 + 1/2) In T + const. 




FIG. 1: (Color online). Form of Aim 2 P„o(t) after time t such 
that z = 2A a t > R 2 , calculated from Eq.flTTJ with z = 2000 
and R — 10 for the three-dimensional walker. The inset for 
lnP n o shows the asymptotic decay. 

Finally, let us note that the strong-field Hamiltonian 
gives similar behaviour. For strong decoherence we 
find that an initially localized state at the origin propa- 
gates as 

/>oo 

Pno(z)= dy e-» P° (zJ M (2Vw)) ■ (14) 
Jo 

and analysis of this shows the same long-time features as 
above. 

Physical interpretation: A path integral analy- 
sis provides some insight here. The anomalous short- 
distance behaviour arises because the effective interac- 
tion between the 2 paths of the density matrix, gener- 
ated by interactions with the spin bath environment, has 
long-time memory effects in it - this is because the bath 
has a degenerate energy spectrum (this is reminiscent of 
weak localisation|2^|). But then how can we explain the 
long-range ballistic tail? Usually even very weak interac- 
tion with a bath gives classical diffusion at long ranges, 
because the environment 'measures' the position of the 
particle as it travels along a given path[23|. For this the 
environment does not have to record all possible trajec- 
tories of the particle - it only needs to track a 'coarse- 
grained' traiectorv|28j. The same is true if the environ- 
ment couples to the particle velocity, from measurements 
of which one can also reconstruct its trajectory. 

The answer to the paradox is interesting. Notice that 
in J3J and JSJ the environmental coupling does not distin- 
guish different particle positions in the space of the graph 
(ie., between different graph nodes), nor the direction of 
transition between them; it only records that transitions 
between them have occurred. This leaves room for the 
constructive interference of many very large paths on the 
graph. 

To gain more insight into the problem we rewrite the 
Hamiltonian i n the momentum representation for the 
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walker and a rotated basis for the TLS spins (rotating 
a k ~* °k)- ^ n this basis the Hamiltonian is diagonal; 
writing <f) = J^k a k< J k, an d given some TLS spin distribu- 
tion {a)*} with a given <j>, then 7i acts on the eigenstates 
|k, {cr|}) according to 

W|k,K}) = cos^ £o (k) |k,K}) , (15) 

If we now start from the initially localized state for the 
walker and arbitrary \\ k |cr|) state for the environment 
(in the original, unrotated basis) the system is an equal- 
weight superposition of all eigenstates. We immediately 
see that states with the same ip evolve coherently with 
a renormalized hopping amplitude A Q cos ip, and in the 
strong coupling limit all values of cp on the [0, 2tt] interval 
are equally represented (thus we rederive the result given 
above in eqtns. JJjJ and 0). The ballistic long-time be- 
haviour comes from those portions of this mixture with 
|cos(^| ~ 1. The anomalous 'sub-diffusive' long-time be- 
havior at the graph origin, on the other hand, comes from 
a small fraction ~ 1/z of states having very small effec- 
tive A cos<p < which cannot propagate anywhere 
at all! 



In quantum information processing systems, where the 
walk can occur in different kinds of information space, 
no general principle forces the environmental couplings 
to distinguish either the different graph nodes, or the di- 
rection of transition between them. Thus we see that in 
the design of quantum computers and certain search al- 
gorithms, it becomes of considerable interest to look at 
quantum walkers for which environmental decoherence 
may even be strong, provided it is not projecting parti- 
cle states onto either the 'position' or 'momentum' bases 
in the information space defined by the graph on which 
the walk takes place. More generally, we see that there 
is an interesting class of systems for which the long-time 
behaviour is very far from diffusive, even in the strong de- 
coherence limit- instead, it combines a short-range 'sub- 
diffusive' behaviour with long-range coherent dynamics. 
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